Influence of measurements on the statistics of work performed on a quantum system 
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The recently demonstrated robustness of fluctuation theorems against measurements [M. Campisi 
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nonequilibrium quantities, such as heat and work, remain unaffected. We determine the impact of 
measurements that are performed during a running force protocol on the characteristic function of 
work. The results are illustrated by means of the Landau-Zener(-Stiickelberg-Majorana) model. In 
the limit of continuous measurements the quantum Zeno effect suppresses any unitary dynamics. 
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I. INTRODUCTION 

The Jarzynski equality and the fluetuation theorems 
are general and surprisingly robust exact relations of non- 
equilibrium thermodynamics. The validity of these rela- 
tions was rigorously established irrespective of the speed 
of the externally applied forcing protocol for driven clas- 
sical [l|-Q and quantum systems 0, Q , staying either in 
complete isolation or in weak Q or even strong contact 
with an environment 0, |1| . See also the reviews 
and references therein. In this work we will focus on the 
Tasaki-Crooks work fluctuation theorem, reading: 
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where pf{w) is the probability density function (pdf) of 
work performed by a forcing protocol denoted as At . This 
forcing acts on the system between the times i = and 
t = T. Accordingly, pb{w) is the pdf of work performed 
on the system when the backward (B) protocol Xr-t de- 
scribes the forcing of the system. The forward and back- 
ward processes start in the Gibbs equilibrium states at 
the same inverse temperature /3 and at the initial pa- 
rameter values Aq and A,-, respectively. The free energy 
difference between these two states is denoted by AF. 

In the derivations 0, 0] of the quantum work fluc- 
tuation theorem, Eq. ([T]), the energy of the system is 
measured at times t = and t = t, and the work w is 
determined by the difference of the obtained eigenvalues. 
Recently, we showed that the Tasaki-Crooks work fluctu- 
ation theorem, Eq. ((!]) , as well as other quantum fluctua- 
tion theorems, remain unaffected for other scenarios than 
this two-measurement scheme: the ratio of forward and 
backward pdf 's in Eq. ([1]) stays unaltered even if further 
projective quantum measurements of any sequence of ar- 
bitrary observables are performed while the protocol is 
in action [ll| . In Sec. |lT]we provide an alternative proof 
of this result, based on the calculation of the character- 
istic function of work of a driven quantum system whose 



dynamics is interrupted by projective quantum measure- 
ments. 

Based on the fact that the process of measurement of 
a quantum system is rather invasive due to the collapse 
of the wave function, it has been argued [ll|, EH that, 
although the value of the ratio of backward and forward 
pdf 's in Eq. ([1} remains unchanged, additional measure- 
ments affect the values of the individual work pdf 's. One 
must expect that, in a driven-measured quantum sys- 
tem, the work done is not only determined by the in- 
teraction with the manipulated external fleld A, but also 
by the measurements themselves, which are physically 
realized by a measurement apparatus. In Sec. IIIII we 
calculate the statistics of work in a prototypical model 
of driven quantum system, namely the Landau- Zener(- 
Stiickelberg-Majorana) [TsHlGj model, and illustrate how 
it is influenced by projective quantum measurements. We 
will also draw the attention on interesting features re- 
lated to the quantum Zeno effect appearing when the 
measurements become very frequent. As we will see, if 
the observable that is measured at time t is the Hamilto- 
nian H{t), then in the Zeno limit the work characteristic 
function approaches the same expression as for an adia- 
batic protocol with no intermediate measurements. 

Sec. IIVI closes the paper with some concluding re- 
marks. 



II. FLUCTUATION THEOREMS FOR 
DRIVEN-MEASURED QUANTUM SYSTEMS 

We consider a quantum system that, in the time span 
[0,t], is thermally isolated, and interacts with the exter- 
nal world only through a mechanical coupling to a time 
dependent external force field and a measurement ap- 
paratus. The information regarding this interaction is 
encoded in a protocol which we denote as 



a = {H{Xt),{U,A)}- 



(2) 
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It specifies a) the system Hamiltonian H{Xt) at each time 
t in terms of the external forces At, and b) the times 
ti G {Q,t), i ~ 1 ... N , a.t which measurements of the ob- 
servables Ai occur. For the sake of simpUcity, in the fol- 
lowing we will adopt the notation H{t) for H{Xt). As in 
the two-measurement scheme 0, 0] of the Tasaki-Crooks 
theorem, wc assume that, besides the N intermediate 
measurements of Ai^ . . . An, the energies determined by 
H{0) and H{t), are measured at times to = and tf — t. 
For simplicity we then set / = + 1 and 



Ao = H{0) Af = H(t) 



(3) 



The remaining observables, Ai,...An, represented by 
hermitean operators are required to have discrete spectra, 
but otherwise can be chosen arbitrarily. 

For times t < the system is assumed to stay in the 
Gibbs thermal state at inverse temperature /3: 



-f3H{a) 



IZO: 



where 



Zo = Tre-^^(O) 



(4) 



(5) 



is the canonical partition function. 

We denote the orthogonal eigenprojectors of the ob- 
servable Ai by n^. These satisfy the eigenvalue equa- 
tions: 



(6) 



where a\ are the eigenvalues of A^i. Hence, with the 
above choice of the first and last measured observables, 
a° and a{ represent the instantaneous eigen-energics of 
the system at t = and t = t, respectively. 

In the following we compute the conditional probabil- 
ity p„{m,T\n,0) to find the eigenvalue obtained in 
the last measurement at time r, provided that the eigen- 
value a° was the result of the first measurement, under 
those conditions that are specified by the protocol a. Wc 
begin our discussion considering only one intermediate 
measurement {N ~ 1) of the observable Ai. 

According to the von Neumann postulates, immedi- 
ately after the eigenvalue a° is measured at ^ = 0, the 
system density matrix becomes: 



p„(o+) = n>X/Pn, 



where 



/Zo 



(7) 



(8) 



is the probability to find the system initially in the state 
with energy a*^. Since wc assumed that the system is 
thermally isolated, it subsequently evolves until time ti 
according to the unitary time evolution Ut-^fl that is gov- 
erned by the Schrodinger equation: 



(9) 



Thus, immediately before the measurement of Ai, occur- 
ring at ii, the density matrix is: 

Pn{tT)=Ut,^OPn{0+)Ulo, (10) 

and the subsequent measurement projects it into 

Pn{t+) =Y,UlUt,,oPn{0+)Ul,nl . (11) 



Likewise, just before the measurement of Af = H[t) at 
time r. the density matrix becomes 



Pn{T ) = Ur.MP7i{tt)Ult^ , 



(12) 



and the probability that the outcome of the measurement 
of Af ~ H{t) at time r is is: 

p„{m,T\n,{)) = TYlilUrMPnittWlt, ■ (13) 

Finally, the pdf of work w, Pa{w), is obtained as the 
sum of the joint probability Pcr(m, T|n, Q)pn restricted to 



(14) 



and hence becomes 

Pc,{w) = '^5[w - a^- Q;°]pcr(m, T|n, 0)p„ . (15) 

A. The characteristic function of work 

Next we focus on the characteristic function of work 
Ga{u), given by the Fourier transform of the work pdf 



G„{u) = / dwP„{w)e' 



(16) 



Substituting into (fT6)) we obtain for the characteris- 
tic function: 

G„{u) = ^e™["--<W(™,r|n,0)e-^</Zo 
= ^ e™["^--"'.lTrnJ:,t/.,t,niC/t„o 

m,n,r 

xulp^^iilul^nlul,^ 



Tr 



e Pq 



(17) 



where [X]^ denotes the time evolution of an operator X 
from t = to t = r in presence of the protocol cr, that 
implies a unitary evolution governed by H(t) interrupted 
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by a measurement of an observable Ai at time t = ti. It 
takes the form 

= E K,o^lultXUr,unlUt,^o, N = l. (18) 



In the case of a protocol a containing N interrupting 
measurements the formal expression of the characteristic 
function is the same as for one interrupting measurement, 
Eq. ([T7|) with the time evolution [X]^ given by 



r/t r/t 



(19) 



B. The Jarzynski equality 



Putting u = one recovers the Jarzynski equality: 



)a = Tr 



Zo' Tr 



Zj/Zo = 



(20) 



because Tr[X]o- — TiX for any trace class operator 
X [13] ■ This follows from the cyclic invariance of the 
trace, the unitarity relation u} JJt^s = 1, and the com- 
pleteness of the projection operators -^r = 1- Here 
Zf = Tre~^^('^) denotes the partition function of the 
Gibbs state at the initial temperature and final parame- 
ter values: 



-liH(T) 



(21) 



and AF = — /3 ^InZo/Zf is the difference of the free 
energies of the thermal equilibrium states pj and Pq. 

The symbol denotes an average with respect to 
the work pdf Paiw). Eq. (PU)) says that the Jarzynski 
equality holds irrespective of the details of the interac- 
tion protocol a. Independent of number and nature of 
the measured observables as well as strength, speed, and 
functional form the driving force the average exponenti- 
ated work is solely determined by the free energy differ- 
ence m. 



C. The work fluctuation theorem 

Besides the Jarzynski equality also the Tasaki-Crooks 
theorem is robust under repeated quantum measure- 
ments. In the presence of many intermediate measure- 
ments it reads: 



■w e 



(22) 



where the tilde (a) indicates the temporal inversion of 
the protocol a, that is: 



{H{T-t),{T-U,A^)}. 



(23) 



Hence, a specifies the succession of force values and 
measurements in the reversed order, specifically with 
the measurement of the observables Ai at times t — ti. 
In particular it implies that at time 0, the observable 
Af ^ H{j) and at time r, the observable = -ff(O) are 
measured. Accordingly, the initial state of the backward 
process is given by the Gibbs state Eq. (PT|) . i.e. the 
system is at equilibrium with inverse temperature /3 and 
force value Ai-. 

Eq. ([22]) holds under the assumptions that both H{t) 
and Ai commute with the quantum mechanical anti- 
unitary time reversal operator [l^. That is, for all 
t e [0, r] , and i = . . . /, we assume 



AS 



(24) 
(25) 



We prove Eq. (1221) for the simplest case of a single 
intermediate measurement, A'^ = 1. The generalization to 
many measurements is straightforward. From Eq. (|17p 
we have 



r 

X nif/i,,oe^^'^-"^^(°V^o. (26) 

Introducing the notation Uf .t for the time evolution gov- 
erned by H{t) = H{t — t), the backward characteristic 
function of work can be written as: 



oe 



--tuH{G) -f}H{T) 



IZf 



=Tr^ e [etc/tt,,o00^nieet[/t ^_,^e 

r 

xett/,_t,,oeete-'("+*«^(")ej eV^/, (27) 

where we used the antiunitarity 09^ = 1. From Eq. 
([25]) it follows that all eigenprojection operators com- 
mute with the time-reversal operator, i.e. 11^8 = 811,^,. 
The time reversal invariance, expressed by Eq. (j24p im- 
plies S^e^'^'-'^O = e/ for any C-number s. Further, 
microreversibility of driven systems [lol . [l^ implies 



Q^Ur-tflQ = Ut,r 



(28) 



Using these relations and recalling that, for any trace 
class operator X, Tr8A8^ = TrA^^, one ends up with: 



G^{u) =Tr^, 



,i{il3-u)H{T) 



X e 



iuH{0) 



KfiKultJZf. 



(29) 
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By comparison with Eq. ((26)) . we finally find 

ZfGa{ip-u) = ZoG,{u), (30) 

hence, by means of an inverse Fourier transform the 
searched fluctuation theorem ((22|) . 

We thus have proved that the fluctuation theorem of 
Tasaki-Crooks remains unchanged if additionally to the 
measurements of energy at time t = and t = t, interme- 
diate measurements of time reversal invariant observables 
Ai are performed at times ti , provided the order of mea- 
surements in the backward protocol is properly changed 
in accordance with the corresponding times t — U. 



III. EXAMPLE 

As mentioned in the introduction, one expects mea- 
surements to strongly influence the pdf of work although 
the Jarzynski equality and the Tasaki-Crooks relation arc 
insensitive to intermediate measurements. To illustrate 
this point in more detail, we consider the example of the 
Landau-Zener(-Stiickelberg-Majorana) p^l - [l6j model de- 
scribed by the Hamiltonian: 

vt 

Hit) = -dz + Aa, . (31) 

It governs the dynamics of a two-level quantum system 
whose energy separation, vt, varies linearly in time, and 
whose states are coupled via the interaction energy A. 
Here, and denote Pauli matrices. 

The Landau-Zcner model is one of the few time- 
dependent quantum mechanical problems that have an 
analytic solution. The elements of the 2x2 unitary 
time evolution matrix Ut^s can be expressed in terms 
of parabolic cylinder functions The instantaneous 

eigenvalues of H{t) arc: 

El^^in- l/2)\/vH^ + 4A2 , n = 0, 1 . (32) 

Since these energies are symmetric with respect to an 
inversion about t = we choose the initial and final times 
to and tf, as ~t/2 and t/2, respectively, instead of and 
T, as in the previous discussion. 

Figure [T^) depicts the survival probability 
p^(0, t/2|0, — t/2) as a function of the instant ti of 
a single intermediate measurement of H{ti) for positive 
times ti and fixed length r of the protocol. As an even 
function of ti this also specifies p^(0, t/2|0, — r/2) for 
negative ti. The straight horizontal line shows the value 
of the survival probability p(0, t/2|0, — t/2) without 
intermediate measurement. It is obvious that the 
intermediate measurement in general alters the survival 
probability thus affecting the work pdf. For example, 
the average work, 

[w) = 2£;[/^tanh[/3£;[/^][l-p(0,T/2|0,-T/2)], (33) 



evidently changes when p(0, t/2|0, — t/2) is replaced by 
p^{0, t/2|0, —t/2). The same can be said for the standard 
deviation of work that reads 

(A«;2) = {2El'y-[l - p{0, t/2|0, -t/2)] - {wf . (34) 

Notably, the introduction of an intermediate measure- 
ment may lower the average work. 

Figure [TJd) shows the N dependence of the survival 
probability p^(0, t/2|0, — t/2) for N equally spaced in- 
termediate measurements of H{t), and two sets of model 
parameters. Oscillatory behavior is observed for small 
values of N, while as N increases the asymptotic value 1 
is approached, see the inset. When the measurement fre- 
quency is high enough the unitary dynamics between sub- 
sequent measurements becomes increasingly suppressed 
until the dynamics is completely frozen and consequently 
the survival probability reaches the asymptotic value 1 



- a) 






— T = n/A.v = 5AVfi. . 



2 4 6 8 10 

ti [A/?i] 




FIG. 1: (Color online) Survival probability of the Landau- 
Zener system, Eq. interrupted by measurements, a) 

Survival probability as a function of the time ti of a sin- 
gle energy measurement of H{ti). The parameter values are 
T = 20ft/A, V — 5A^/ft. The horizontal line is the survival 
probability p(0, — r/2|0, r/2) with no intermediate measure- 
ments, b) Survival probability as a function of the number 
A*' of equally spaced intermediate measurements of H{t). The 
parameter values are r — 20h/A, v = 5A^/h (black circles), 
and T = 10ft/ A, v = 20A^/ft (red triangles). Inset: the sur- 
vival probabilities approach the value 1 in the Zeno-limit. 
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for iV ^ oo. This phenomenon is known as the quan- 
tum Zeno effect |2l|-|24|. We investigate it further in the 
foUowing section. 



A. Quantum Zeno Effect 

To formaUy elucidate the quantum Zeno effect ob- 
served in this particular example, and under more general 
conditions as well, we analyze the form of the character- 
istic function given by Eq. (jl7p in the limit of infinitely 
many measurements of energy. We approach this limit 
by considering a finite number of N intermediate mea- 
surements of energy that take place at equally spaced 
instants 



tk=to + ke , 



k = l...iV, 



(35) 



where e — {tf — to)/{N + 1) denotes the time elapsing 
between two subsequent measurements. We denote the 
corresponding protocol with symbol v, that is: 



iy = {H{t),{tk,H{tk))}. 



(36) 



Since we are interested in the limiting case of infinitely 
many measurements, we may choose N sufBcicntly large 
such that the Hamiltonian between two subsequent mea- 
surements can safely be approximated by its value at the 
later measurement, i.e. H{t) w H{th) for t e (tfc_i,ifc) 
with tfc = to + ke being the instant of the fcth measure- 
ment. The time evolution within such a short period then 
becomes 



t 



(37) 



acts 



According to Eq. ([T9l) the time evolution 
on the projection operator 11^^. (tfc) resulting in the phase- 
factor 6*"'"'= while the complex conjugate phase fac- 
tor is obtained from the product Iirk(tk)Ut 



tk.tk- 



which 



appears right of the operator X [30[ • Hence, these factors 
cancel each other and one obtains for the :^-propagated 



exponential operator X 



AuHir) 



the expression 



J2 TlrAtl)Tlr,{h) . 



ri,r2....rN 



xn,„(tjv)e™''^"^n,„(tAr)., 



(38) 



Assuming first that H{t) is non-degenerate at any time 
t G (tQ,tf) all projection operators are of the one- 
dimensional form Ilr^{tk) = |'0rfc(tfc))(V'rA.(tfc)| with in- 
stantaneous cigcnfunctions \iprkitk)) of the Hamiltonian 
H{tk) ■ The labeling of the instantaneous eigenstates can 
be arranged in such a way that the eigenvalues ak (t) are 
continuous functions of time. In other words, each adia- 
batic energy branch is labeled by an index r. The scalar 
products of cigcnfunctions at neighboring times then de- 
viate from a Kroncckcr delta by terms of the order 



which consequently can be neglected, i.e., states on dif- 
ferent adiabatic energy branches at neighboring times are 
almost orthogonal 



(39) 



Herewith the left hand side of Eq. (j38|) simplifies to read 
in the limit of infinitely many dense measurements 



E 



\Mto)){MtfW"^''^\Mtf)){Mto)\- 

(40) 



For the generating function this yields 



,m[a,(tj)-Q,(to)]„-/3a,(to) 



/Zo . (41) 



We therefore find that in the Zcno-limit the characteristic 
function of work and accordingly the pdf of work coincide 
with the respective expressions obtained for an adiabatic 
protocol in the absence of any intermediate measurement. 

The same line of reasoning also applies for the case that 
an observable that does not change in time is repeatedly 
measured, i.e., for a protocol 



^i = {H{t),{tk,A)}. 



(42) 



For such an observable A with eigenprojection operators 
Ilr one then obtains in the Zeno limit for the character- 
istic function 



G'^^iu) = Tr^ (Ure'""^''^U, 



g-mH(to)g-/JH(to) 



(43) 



1. Landau- Zener 



In the case of the Landau-Zener problem studied in 
the previous section, Eq. (|3ip . the eigenvalues on each 
adiabatic branch coincide at the beginning and the end 
of the protocol and therefore the coefficients of u in the 
exponential terms on the right hand side of Eq. (PT|) 
vanish leading to 



1 



(44) 



This leads to the expected result that under permanent 
observation of the Hamiltonian no work is done, i.e.. 



(45) 



The Zeno effect sets in when the measurement frequency 
is so large that any unitary evolution within the time 
leaps between two subsequent measurements can be ne- 
glected, hence for e^Jv^t"^ + 4A/fi, ^ 1, or cquiva- 
lently if TV > tsJv^t"^ + 4A/?i. With the parameters 
used in Fig. [TJj), the Zeno effect sets in for iV ^ 10'^, as 
shown in the inset of Fig. [Ud) . 

In the case of degeneracy. Eq. (PT|) continues to hold if 
the instantaneous energy branches arc labeled such that 
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the corresponding instantaneous eigenvalues smoothly 
vary along each branch. For the Landau-Zener problem 
a degeneracy happens at t = if the coupling strength 
A vanishes. One then finds the characteristic function 

cosh(pur/2) 

yielding for the pdf of work 

P;^{w) ^p^S{w~vT)+p+d{w + vT) [A = 0], (47) 

where p± = e^F/^^^/^/p cosh(/3T/2)] denote the thermal 
populations of the ground- and the excited states at the 
initial time t ~ —t/2, respectively. 

In the case of countinuous measurement of the Pauli 
operators ct;, i = x,y,z, the driving-measuremct proto- 
cols are 

= (t/c,(7i)}, i^x,y,z. (48) 

The characteristic functions of work then become 

1 / A2 



(49) 



G.-H = i(i + f)(.- 



where p± denotes the populations of the exited and the 
ground states of the initial Hamiltonian, respectively, 
given by 



with 



v/(W2)2 + A2 



(50) 



(51) 



The work pdf follows then by means of an inverse Fourier 
transform to read 



pZ ^ ^iP-^i^ - 2?) +P+6{^ + 2g) + 5{w)) (52) 
A2\ 

1 + — {p^5{w - 2q) +p+d{w + 2q)) 



1-^1 5{w) . 



IV. DISCUSSION 

We provided an alternative proof of the robustness of 
the Tasaki-Crooks fluctuation theorem against repeated 
measurements The present proof however is more 

general in regard to the fact that it allows for possi- 
bly degenerate eigenvalues of the measured observables. 
Further it was found that the characteristic function of 
work assumes the form of a quantum two-time correlation 
function of the exponentiated initial Hamiltonian, and 
the exponentiated final Hamiltonian in the Heisenberg 
picture generated by the unitary evolution interrupted 
by quantum collapses. 

Our proof keeps holding also if quantum measurements 
are done with respect to positive operator valued mea- 
sures (POVM) [25 which are less invasive than those 
described by von Neumann projection valued measures. 
For the case of measurements with respect to a POVM, 



^ MrM}. = 1 



(53) 



the interrupted evolution \X\cj is obtained from Eq. (fT9|) 
with the projection operators H^ on the left hand side 
of X being replaced by the weak measurement operators 
Mr^ and those right of X by the adjoint operators Aft. 
Due to Eq. ([521), Tr[X]^ = TrX continues to hold imply- 
ing the validity of the Tasaki-Crooks relation, Eq. (|22p . 
for weak measurement of time reversal invariant observ- 
ables with MfeO = OMfe. 

In proving Eq. ([22|) we assumed that the Hamiltonian 
and the measured observables commute with the time 
reversal operator. This assumption may be relaxed and 
the Tasaki-Crooks fluctuation theorem continues to hold 
if the backward protocol is defined as: 



b = {ei/(T - t)Q\ (t - t„ QA^)} 



(54) 



Notwithstanding the robustness of the Tasaki-Crooks 
theorem, we found, in accordance with the intuitive ex- 
pectation, that the interruption of the dynamics of a 
driven quantum system by means of projective measure- 
ments alters the statistics of work performed on the sys- 
tem. We illustrated the infiuence of measurements on the 
work distribution by the example of the Landau-Zener(- 
Stiickelberg-Majorana) model. We noticed that depend- 
ing on the positions in time or frequency of one or more 
interrupting measurements, the average work {w) may be 
lo wered. Fig. [H 

ISchmiedl et al.l [26j studied the problem of design- 
ing optimal protocols that minimize the average work 
spent during the forcing, without considering interme- 
diate measurements. In order to further minimize the 
average work, one might expand the optimization pa- 
rameter space and include the possibility of performing 
intermediate measurements. This opens for the possibil- 
ity of a much greater control over the energy flow into 
the system. 

In the limit of very frequent measurements the uni- 
tary dynamics becomes completely suppressed due to 
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the quantum Zcno effect. If at each, measurement time 
the instantaneous Hamiltonian is measured, then in the 
Zeno hmit the work characteristic function approaches 
the same form that it would assume if the protocol was 
adiabatic. 

The result of an intermediate measurement may be 
used to alter the subsequent force protocol. In this way 
a feedback control can be implemented for classical [27 1 
as well as for quantum systems In both cases the 

Jarzynski equality only holds in a modified form. 

Landau- Zener dynamics and frequent quantum mea- 
surements of a two-level system coupled to a thermal 
bath were recently shown to provide efficient means for 
quantum-state preparation and purification [2§|, re- 



spectively. These are crucial prerequisite for the imple- 
mentation of working quantum computers. The com- 
bined study of dissipative Landau-Zener dynamics with 
frequent observations, could unveil yet new practical 
methods for quantum state control and manipulation. 
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